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APPENDIX  I 


F-W  EXTENDED  ANQLE-OF-AITACK 
SIMULATION  MODEL 


A.  INTRODUCTION 

This  appendix  documents  the  mathematical  model  used  in  the  unpiloted 
(digital)  and  piloted  high  AOA  simulations  of  a  generic  high-speed,  twin-jet, 
fighter- type  airplane.  The  model  represents  a  relatively  "clean"  aircraft 
(no  flap,  slat,  or  gear  deflection)  based  upon  an  F-4j,  with  external  stores 
consisting  of  one  centerline  fuel  tank  and  one  air-to-air  missile .  Four 
different  generic  configurations  are  generated  by  changing  a  limited  number 
of  aerodynamic  coefficients. 

The  appendix  is  organized  as  follows:  Section  B  presents  the  equations 
of  motion  including  a  discussion  of  assumptions,  axis  systems  and  transfor¬ 
mations.  The  equations  for  the  aerodynamic  force  and  moment  coefficients 
are  presented  in  Bection  C.  The  airplane  physical  characteristics  and  con¬ 
trol  system  are  described  in  Section  D.  Section  E  summarized  the  data  sources. 
The  aerodynamic  coefficient  data  are  presented  in  Section  F.  The  coefficient 
"lookup"  tables  are  listed  sub sequentially  in  Table  9  and  plotted  in  Figs.  5 
to  20.  Validation  of  the  model  with  flight  test  data  is  given  in  Section  0. 

B.  EQUATIONS  07  MOTION 

The  nonlinear  slx-degree-of-freedom  aircraft  equations  of  motion  used  in 
the  simulation  programs  are  consistent  with  the  following  assumptions. 

1 .  Assumptions 

a.  The  airframe  is  assumed  to  be  a  rigid  body. 

b.  The  earth  is  assumed  to  be  fixed  in  inertial  space. 

c.  The  mass  and  mass  distribution  of  the  vehicle  are 
assumed  to  be  constant. 

d.  The  aircraft  has  a  plane  of  symmetry. 
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e.  Effects  associated  with  rotation  of  the  vertical 
relative  to  inertial  space  are  assumed  negligible; 
the  magnitude  of  the  gravity  vector  is  assumed 
constant. 

2.  Axis  Systems 

All  the  axis  systems  (see  Pig.  i)  have  their  origin  at  the  aircraft 

center  of  gravity  (c.g.)«  The  Earth  Axis  system  is  oriented  with  the  ZE 
axis  along  the  local  gravity  vector  and  the  Xg,  axes  in  a  horizontal 
plane  with  arbitrary,  but  fixed,  direction. 

The  Body  Axes  are  fixed  in  the  aircraft  with  the  Xg  axis  positive 
forward  along  the  fuselage  reference  line.  The  Zg  axis  is  in  the  plane 
of  symmetry,  positive  down,  and  the  Yg  axis  is  perpendicular  to  the  plane 
of  symmetry,  positive  out  the  right  wing.  The  Body  Axes  are  located 
relative  to  the  Earth  Axes  by  conventional  Euler  Angles  Y,  0,  and  i>. 

V  is  a  rotation  of  the  Body  Axes  from  the  Earth  Axes  about  the  Zq  axis; 

8  is  a  rotation  about  the  intermediate  Yg  axis;  and,  finally,  <t>  is  about 
the  Xg  axis.  All  angles  and  angular  rates  are  positive  in  a  righthanded 
sense. 

The  final  axis  system  is  called  Flight  Path  Axes.  They  are  aligned 
with  the  aircraft  inertial  velocity.  For  the  still-air  case  considered 
here  they  are  identical  to  Wind  Axes  (usually  aligned  with  the  aircraft 
velocity  relative  to  the  air  mass).  Herein  we  use  the  more  conventional 
W  subscript  and  Flight  Path,  Flight,  and  Wind  interchangeably.  The  axis 
lies  along  the  aircraft  total  velocity  vector,  Vj;  the  Zyr  axis  is  in  the 
aircraft  plane  of  symmetry;  and  the  Y ^  axis  completes  the  right-handed 
orthogonal  set.  The  Wind  or  Flight  Axes  are  located  relative  to  the  body 
by  the  angles  of  attack,  a,  and  sideslip,  £;  a  is  a  rotation  about  the 
axis  and  £  a  rotation  about  the  Zw  axi-3  (if  a  =  0) . 

The  transformation  from  Flight  Path  Axes  to  Body  Axes  is  given  by: 


V 

cos  a  cos  £ 

-cos  a  sin  £ 

-sin  a 

XW 

a 

sin  £ 

cos  £ 

0 

YW 

ZB 

sin  a  cos  £ 

-sin  a  sin  £ 

cos  a 

a 

. 

_ 
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3 ,  Moment  Equations 


P>  Q,  and  R  are  the  instantaneous  components  of  the  aircraft  total 
inertial  angular  velocity  vector  in  the  Xjj,  Y$j,  Zg  Body  Axes.  The  deriva¬ 
tives  of  these  components  with  respect  to  time  are  related  to  the  applied. 

moments  about  the  center  of  mass  by: 


where 

C1 

c2 

c3 

c4 

°5 


• 

P  - 

(°1®  + 

CgP)Q 

+  c3  X, 

Q  - 

e^RP  + 

C6^r2 

-  Ps) 

+  Cj&L 

(2) 

R  - 

(ceP  + 

CgR)^ 

+ 

+  c10& 

°{(v 

-  Iz)^z 

-  iiz} 

eg  » 

Ixz/Iy 

°{  V 

+  : 

Ixz 

Qij  m 

1  ./Iy 

(3) 

0  iz 

00  «= 

o{  (Ix-  Iy)lx  +  4Z) 

G  Ixz 

0<j  « 

G  {  Iy  ”  ”  I*  |  Ixz 

H 

1 

N 

M 

x)/!y 

c10  “ 

OIx 

a  -  i/dxiz  -  ix2 ) 


Ix,  ly,  Iz,  and  Ixz  are  the  moments  and  product  of  inertia  with  respect  to 
the  Body  Axes.  £ ,  YYl,  and  7L  are  moments  about  the  Body  Axes  due  to 
aerodynamics  and  aircraft  thrust. 


4.  Force  Equations 

Flight  Path  Axes  force  equations  are  used  for  computational  efficiency. 
The  derivatives  of  a,  p,  and  the  aircraft  total  inertial  translational  veloc 
ity,  VT,  are  related  to  the  external  forces  by: 


ati  -va'i*  IK'.iV  .  \,f!  -i,\ M.tt  l/cJibutLtlx 


a  «*  Q,  —  tan  p(P  cos  a  +  P.  sin  a)  +  Z^/(mVip  cos  P) 

P  =  P  sin  a  -  R  cos  a  +  Yw/(mVT)  (Vi 

Vj  a  3fy/m 

XW,  Y^,  and  Zy  are  the  components  of  the  total  external  forces  (aerodynamic, 
thrust,  and  gravity)  along  the  Flight  Path  Axes*  m  is  the  aircraft  mass* 

The  a  equation  is  limited  in  accuracy  to  small  sideslip  angles.  The 
large  excursions  in  angles  of  attack  and  sideslip  encountered  in  the  piloted 
simulation  necessitated  a  change  in  definition  of  a*  The  Body-Axis  equation 
for  vertical  acceleration,  w,  was  substituted: 

V'  *=  QVrp  cos  p  -  Vrji  sin  p  (P  cos  a  +  R  sin  a)  +  Z^/m 

and  the  definitions  for  a  and  0  were  a  «  sin“1  w/Vrp,  p  =  sin“"'  v/Vj* 
These  excursions  were  not  experienced  in  the  preliminary  analytical  work  per 
formed  in  preparation  for  the  piloted  simulation* 

3 .  Thrust  Geometry 


The  thrust  is  assumed  to  lie  in  the  aircraft  plane  of  symmetry  and  is 
oriented  with  respect  to  the  Body  Axes  by  the  thruet  inclination,  if  and 
offset,  Zj,  as  indicated  in  the  sketch  below. 


x» 


6,  Body  Axis  Moments 

The  roll, ,2^  ,  and  yaw,?L,  moments  are  due  to  aerodynamics  only,  while 
the  pitching  moment,  771,  includes  a  thrust  term.  These  are  given  by; 


oC  =  qSbC^  f 


qScCm  +  ZjT 


)  %  »  qSbCn 


(5) 


. - . 


is  the  dynamic  pressure;  b,  £,  and  S  are  the  reference  span,  chord,  and 
wing  area,  respectively.  C£,  Cm,  and  Cn  are  the  total  body-axis  non- 
dimensional  aerodynamic  moment  coefficients  (referenced  to  the  c.m.). 

7.  Flight  Path  Axis  Forces 


The  forces  along  the  Flight  Path  Axes  have  components  due  to  thrust  (T), 

weight  (mg),  and  aerodynamics.  They  are  given  by: 

V 


T  cos  p  cos  (a  +  |0) 

+  mgjcos  0  cos  4>  sin  a  cos  P  -  sin  ©  cos  a  cos  P  +  cos  0  sin  $  sin  P 
+  Ya  sin  p  -  D  cos  p 


-T  sin  P  cos  ( a  +  40) 

+  mgjcos  0  sin  $  cos  P  +  sin  0  cos  a  sin  p 
+  YA  cos  p  +  D  sin  p 


-  cos  9  cos  <t  sin  a  sin  p 


Zy  «  -T  sin  (a  +  £Q)  +  mgisin  0  sin  a  +  cos  0  cos  ®  cos  a>  -  l 


L,  D,  and  YA  are  the  aerodynamic  forces  and  are  given  by: 

L  *  3SCL 

d  =  -qscD  (7) 

YA  -  gSCy 

CL,  CD,  and  Cy  are  the  lift,  drag,  and  side  force  coefficients.  Lift,  L, 

and  drag,  D,  are  assumed  in  the  plane  of  symmetry;  lift  positive  up  and  per¬ 
pendicular  to  the  velocity  vector,  drag  positive  aft  and  along  the  projec¬ 
tion  of  the  total  velocity  in  the  plane  of  symmetry.  The  side  force  YA  la 
that  component  of  the  total  aerodynamic  force  which  is  perpendicular  to  the 
plane  of  symmetry;  it  lies  along  the  YB  axis  and  is  positive  out  the  right 
wing. 
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It  should  be  noted  that  the  drag  and  side  force  are  not  the  same  as  the 
X  and  Y  Flight  Path  Axes  forces  for  non-zero  sideslip,  see  Fig.  2.  The  aero¬ 
dynamic  forces  instead  have  been  defined  in  terms  of  the  axes  in  which  aero¬ 
dynamic  data  are  most  frequently  available. 

8 i  Body  Axis  Euler  Angle  Rates 

The  Body  Axis  Euler  angle  rates  are  related  to  the  angular  velocity 
components  by: 


T  ■  (R  cos  0  +  Q  sin  0)/cos  © 

8  ■  Q  eos  0  -  R  sin  0  (8) 

9  ■  P  +  f  sin  8 

9.  Body  Axis  Velocities 

The  components  of  the  total  velocity  along  the  Body  Axes  are: 

U  ■  Vrp  cos  a  cos  p 

V  -  VT  sin  P  (9) 

W  ■  VT  sin  a  cos  P 

10.  Earth  Axis  Velocities 

The  Body  Axis  velocities  are  transformed  through  the  Euler  angles  to 
yield  the  Earth  Axis  velocities. 

Xg  ■  U  cos  0  cos  ▼  +  W(eos  0 

+  V(sin  <t 

ye  ■  U  cos  0  sin  y  +  v(sin  9 

+  W(cos  <t> 

Zg  ■  -U  sin  0  +  V  sin  9  cos  0  +  W  cos  9  eos  0 
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sin  0  cos  v  +  sin  9  sin  y) 

sin  0  cos  v  —  cos  0  sin  y) 

sin  0  sin  V  +  cos  9  cos  v)  (10) 

sin  0  sin  y  —  sin  9  cos  v) 


'  -  J.iSh  UlV  yt  L'_L>  !*•_ .  1 
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Figure  2.  Relation  of  Aerodynamic  Force  Axis  Syst 

to  Wind  Axes 


11.  Aircraft  Controls 

The  aircraft  is  assumed  to  have  pitch,  roll,  and  yaw  controls,  i.e., 
®stab>  ®a,  ^sp>  ant*  ^r*  Positive  deflections  and  travel  limits  for  these 
controls  are  given  in  Section  D  of  this  appendix. 

12.  Auxiliary  Variables 

Additional  variables  which  are  useful  in  aircraft  dynamics  studies 
are  given  below. 

a.  Flight  Path  Inclination  Angle,  v 


tan  7  -  “Z/n/x^  +  Y2 
sin  7  “  -Z/Vj 


sin  7  *  sin  8  cos  a  cos  3 
-cos  8  sin  4  sin  3 
-cos  8  cos  4  sin  a  cos  3 


b.  Altitude  Rate,  K 


H  -  -Z 


c.  Body-Mounted  Accelerometer  Signals,  at  Center  of  Gravity 
axcg>  aycg>  azcg 

r  ax  T  Uo«  a  cos  3  —cos  a  sin  3  -sin  a  "1 


sin  3 


cos  3 


sin  a  cos  3  -sin  a  sin  3  cos  a 


X  Vip(3  -  P  sin  a  +  R  cos  a) 

Vj  cos  3[&  ■  Q  *•  tan  3(P  cos  a  +  R  sin  a)] 


-sin  8 


-  g|  cos  6  sin  4 


cos  6  cos  4 


. y.n;, (Si,. V.  VI in. 


or 


T  cos  |0  —  D  cos  a  +  L  sin  a 
*A 

~T  gin  S0  “  D  sin  a  ”  L  coa  a 


d.  Body-Mounted  Accelerometer  Signals  at  Arbitrary  Point 
(ix,  Jy#  ia)  in  Body,  a^,  *Jr,  a£ 


Zx,  and  are  the  coordinates  of  the  arbitrary  point  in  Body  Axes, 
e.  Stability -Axis  Angular  Velocities,  «s>  *a 

A  body-fixed  Stability  Axes  system  (Ref.  l)  may  be  defined  which  is 
located  with  respect  to  Body  Axes  by  the  trim  angle  of  attack,  aQ.  The 
components  of  the  aircraft's  total  inertial  angular  velocity  in  body-fixed 
Stability  Axes  are  given  by: 

Pa  =>  P  cos  a0  +  R  sin  a0 

Qs  B  ft  (15) 

Rs  «  -P  sin  a0  +  R  cos  aQ 


f.  St  ability- Axis  Euler  Angles,  Yg,  ©g,  ©s 

The  body-fixed  Stability  Axes  system  defined  above  can  be  located  relative 
to  Earth  Axes  by  a  conventional  Euler  Angle  set,  Ys,  ©s,  ©g.  The  Stability 
Axis  Euler  Angles  are  related  to  the  Body  Axes  Euler  Angles  and  the  trim 
angle  of  attack  by: 


tan  Y. 


cos  a0  sin  Y  cos  0  +  sin  ao(sin  Y  sin  0  cos  ©  -  cos  Y  sin  ©) 
cos  a0  cos  Y  cos  0  +  sin  a^cos  Y  sin  ©  cos  $  +  sin  Y  sin  ©) 


sin  ea  =  cos  a0  sin  0  -  sin  a-  cos  0  cos  0 


tan  0g  -  cos  ©g 


cos  a©  sin  0  —  sin  a0  cos  0  cos  © 
cos  a0  cos  0  cos  ©  +  sin  a0  sin  0 


cos  9  sin  4 


tan  © 


s  cos  a©  cos  0cos  ©  +  sin  Oq  sin 


g.  Stability-Axis  Euler  Angle  Rates,  Ys,  0a,  ®E 


ts  “  (Ra  C0S  ®s  +  ^s  ain  ^s)/cos  ®s 
ea  ■  Qg  cos  ©a  -  Rg  sin  ©g 
©g  B  Ps  +  Yg  sin  0g 

h.  Flight  Path  Azimuth  Angle,  X 

tan  1  =  Y/X 


C.  AERODYNAMIC  FORCE  AND  MOMENT  COEFFICIENT  EQUATIONS 


Aerodynamic  forces  and  moments  have  been  defined  in  terms  of  the  total 
non-dimensional  aerodynamic  coefficients,  i.o.,  the  stability-axi*.  force 
coefficients  C^,  C^,  Cy  and  the  body-axis  moment  coefficients  0^,  Cm,  Cn. 

As  shown  in  Fig.  5,  these  are  functions  of  the  aircraft  state,  the  control 
surfaces,  tabular  coefficient  and  derivative  data,  and  constant  parameters. 
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Coefficients  and  Derivatives 


The  coefficient  equations  are  given  in  Table  1 .  A  summary  of  the  tabu¬ 
lar  functions  is  given  in  Table  2  (the  data  is  given  in  Section  F)  and  the 
constant  parameters  are  presented  in  Tables  3  and  5.  As  indicated  in  Table  2, 
all  tabular  data  are  given  for  p-deg  increments  in  a  for  either  0  deg  <  a  < 

110  deg  or  0  deg  <  a _<  45  deg.  Linear  interpolation  between  adjacent  data 
points  was  used  to  obtain  derivative  values  for  intermediate  values  of  a. 

All  coefficients  and  derivatives  except  ClbasiC  ahd  ^BASIC  8x6  symmetric 
functions  of  a  and  the  absolute  value  of  a  was  used  to  obtain  the  derivative 
value.  Cm  is  an  anti- symmetric  function  of  a;  as  indicated  in  Table  1, 

|a|  is  used  to  obtain  the  table  value  and  the  sign  of  the  appropriate  term 
is  changed  for  negative  a' s.  In  general  for  |a|'s  which  exceed  the  range 
of  tne  tables,  the  derivative  value  for  the  maximum  |a|  given  in  the  table 
is  used.  This  is  shown  explicitly  in  Table  1,  where  the  complete  equations 
have  been  given  for  each  a  range.  It  should  be  noted  that,  except  for  Cl, 
the  equations  themselves  are  not  dependent  on  the  a  range.  is  not 

quite  anti- symmetric;  the  explicit  equations  for  lift  coefficient  at  nega¬ 
tive  angles  of  attack  are  given  in  Table  1 . 

ACmit  is  the  only  two-dimensional  function  used;  the  table  contains  values 
for  |  p|  from  0  deg  to  30  deg  in  5  deg  increments.  For  |p|  >  30  deg,  the  table 
values  for  p  >=  'IQ  deg  were  used.  The  organization  of  the  data  1b  discussed 
further  in  Section  F. 

As  Table  2  indicates,  there  are  two  sets  of  data  for  each  of  bhe  deri¬ 
vatives  c£pj>  a^d  Selecting  the  proper  combinations  of 

these  derivatives  results  in  the  desired  four  basic  aircraft  configurations 
as  indicated  in  Table  b- . 
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TABLE  1  .  AERODYNAMIC  COEFFICIET5T  EQUATIONS 
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TABLE  1.  AERODYNAMIC  COEFFICIENT  EQUATIONS  (Concluded) 


b)  Body  Axis  Moment  Coefficients 
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TABLE  2.  ORGANISATION  OF 
TABULAR  AERODYNAMIC  DATA 


PARAMETERS 

FUNCTION 

OF 

UNCT3 

MNEMONIC 

AND  ARRAY 
DIMENSIONS 

LINE 

NUMBERS 

a 

ORID# 

Cl3ASIC 

a 

- 

CLBAS( 23) 

2-8 

o.(5.)iio. 

ClSstae 

a 

l/dag 

CLSTAB(23) 

9-15 

o.(5.)lio. 

CdBASXC 

a 

- 

CDBAS(23) 

16-22 

0.(5.  )no. 

% 

a 

l/deg 

CYB(SJ) 

23-29 

0.(5. >1 10. 

°yBr 

a 

l/deg 

CYDR(lO) 

30-33 

0.(5. )43. 

°-p, 

a 

l/d*g 

CRB1 ( 23) 

34-40 

0.(5. )110. 

c4p2 

a 

l/dag 

CRB2('23) 

41-47 

0.(5.)U0. 

c'p, 

a 

l/rad 

CRPI(IO) 

48-51 

o.(5.)45. 

cipg 

a 

l/rad 

CRP2(10) 

52-55 

o.(5. )45. 

°<r 

a 

l/rad 

CRR(lO) 

56-59 

0.(5045* 

°<». 

a 

l/dag 

CRDA(lQ) 

6o-63 

0.(5.)45. 

C4tjgp 

a 

l/dag 

CRDSP(lO) 

64-67 

o.(S.)45. 

C/Br 

a 

l/dag 

crdh(io) 

68-71 

0.(5. )45. 

C«BASIO 

a 

- 

CMBAS(23) 

72-78 

0.(3. )  1 1 0. 

ACmj 

a,  0 

- 

DCMl(lO,7) 

79-95 

0.(5. )45. 

A0m2 

a,  P 

■M 

DCM2(l0,7> 

96-112 

0.(5.)45. 

°mq 

a 

l/r«d 

OMQ( 1 0) 

113-116 

0.(5. )45. 

am& 

a 

l/rad 

cmad(io) 

117-120 

o.(5- >45. 

CmtSTAB 

a 

l/deg 

CMSTAB(25) 

121-127 

o.(5.)lio. 

C|% 

CL 

l/deg 

cmda(io) 

128-131 

0.(5.145. 

Cm6ap 

a 

l/deg 

cmdsp!io) 

132-135 

0.(5045. 

CnPl 

a 

1/deg 

CNB1  ( 1 0) 

136-139 

0.(5. )45. 

CnP2 

a 

l/deg 

CNB2(l0) 

l4o-i43 

0.(5. )45. 

cnp 

a 

l/rad 

CNP{  10) 

1 44-1 47 

0.(3-)45. 

Cnr 

a, 

l/rad 

CNH(lO) 

148-151 

o.(5-)45. 

C% 

a 

l/deg 

cnda(io) 

152-155 

0.(5.)45. 

Cru 

°op 

a 

l/deg 

CNl)SP(lO) 

156-139 

0.(5. )4j. 

% 

°r 

a 

l/deg 

cram(io) 

1 60-1 63 

0.(5.145. 

*  0.(9.)110.  indicated  valued  ore  given  for  a  from  0  dog  t,o  110  dog 
In  5  dug  Increment*,  l.o.,  0°.,  5°.,  10°,,  .  lot)0.,  109°, ,  110°. 
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TABLE  3.  CONSTANT  AERODYNAMIC  COEFFICIENTS 


ACdstores^ 

.0037 

Cy6Jl/deg) 

-.0001 67 

Cy_  (l/deg) 

Bsp 

-.00006 

^Bstab(1/deg2) 

.00011 

TABLE  4.  AIRCRAFT  CONFIGURATTON/DERIVATIVE  COMBINATIONS 


AIRCRAFT 

DERIVATIVE  INDEX 

CONFIGURATION 

i 

\ 

J 

ACmk 

k 

CnP£ 

i 

A 

1 

1 

1 

1 

B 

1 

2 

1 

1 

C 

2 

1 

1 

1 

D 

1 

1 

2 

2 

S.  AIRPLANE  PHYSICAL  CHARACTERISTICS  AMD 
CONTROL  SYSTEM 

1 .  Physical  Characteristics 

This  section  describes  the  pertinent  dimensions  of  the  airplane  (weight, 
wing  area,  inertias,  etc.)  and  defines  the  control  surfaces.  The  model  is 
based  upon  the  F-4J,  but  is  considered  to  be  representative  of  a  generic 
fighter-type  airplane . 

The  major  dimensional  parameters  for  this  model  are  listed  in  Table  5, 
the  control  surfaces  are  outlined  below  and  control  surface  limits  are  given 
in  Table  6. 

2.  Control  System 

Longitudinal  control  is  provided  by  an  all-moving  horizontal  tail.  Lat¬ 
eral  control  is  provided  by  a  combination  of  spoilers  and  ailerons.  The  ail¬ 
erons  deflect  downward  only;  the  spoilers  deflect  upward  only.  The  left  ail¬ 
eron  and  right  spoiler  operate  simultaneously,  as  do  the  right  aileron  and 
left  spoiler.  Spoiler  and  aileron  deflection  are  combined  by  a  simple  relations 

&sp  "  1*433&a 

Directional  control  is  provided  by  a  conventional  rudder. 

The  longitudinal  control  system  is  outlined  in  Table  7,  and  the  lateral/ 
directional  control  system  is  shown  in  Fig.  4. 


TABLE  5.  AIRPLANE  DIMENSIONAL  PARAMETERS 


SYMBOL 

DEFINITION 

VALUE 

Id 

Wing  span,  ft 

58.67 

C 

Wing  mean  aerodynamic  chord,  ft 

16. 04 

S 

Wing  area,  ft® 

550 

w 

Airplane  weight,  lb 

37,000 

xref 

Reference  center  of  gravity, 

31 

xcg 

Actual  center  of  gravity,  #c 

29.3 

Ix 

Moment  of  inertia  about  X-body  axis, 

slug-ft2 

23,850. 

Iy 

Moment  of  inertia  about  Y-body  axis, 

slug-ft2 

127,^00. 

Moment  of  inertia  about  Z-body  axis, 

slug-ft2 

1)4.6,000. 

Ixz 

Product  of  inertia  about  X-Z  axes, 

slug-ft2 

2210- 

Thrust  inclination,  dcig 

5*25 

zd 

Thrust  offset,  ft 

-0.336 

H 

Reference  altitude,  ft 

15,000 

19 


TABLE  6. 

CONTROL  SURFACE  DEFINITIONS 


SURFACE  DEFLECTION  LIMITS ,  DEG 


Bsta'b>  Positive  TED' 

Ba,  Positive  Left  TED 
Bap,  Positive  Right  TEU 
6r,  Positive  TEL 


21  up,  9  clown 
0  up,  30  down 
43  up,  0  down 
430 


TABLE  7. 

LONGITUDINAL  CONTROL  SYSTEM 


•  No  SAS  or  CAS 

•  Basic  F-4  manual  control  system  gain, 

Bstab  slong  =  dee/in- 

•  Feel  system  —  simple  spring/damper 

(gradient  *^3.3  lb/in) 


Breakout  —  1 .5  to  2  lb 
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S.  DATA  SOURCES 


This  section  outlines  the  sources  for  the  aerodynamic  data  presented  in 
this  appendix.  Reference  1 6  was  used  as  the  initial  data  base  for  this  simu- 
lation  model.  The  Ref.  16  data  covers  an  a  range  of  -5  deg  to  31  deg  and  a  p 
range  of  —10  deg  to  10  deg.  Wind  tunnel  tests  of  an  F-4  model  (Ref.  2),  cov¬ 
ering  an  a  range  of  -10  deg  to  110  deg  and  a  p  range  of  -40  deg  to  40  deg., 
were  the  primary  sources  for  the  extension  of  static  coefficients.  Dynamic 
coefficients  were  extrapolated  based  on  trends  of  data  available  from  other 
sources  to  60  deg  angle-of-attack,  with  the  intent  of  making  such  coefficients 
simple  to  mechanise . 

Because  of  the  nature  of  the  reference  (an  STI  working  paper),  some  dis¬ 
cussion  is  in  order  on  the  data  which  was  used  to  develop  the  initial  data 
set.  Three  separate  sources  were  utilized  in  Ref.  16. 

The  first  two  sets  of  data  were  supplied  to  STI  for  a  previous  study. 

One,  Ref,  3,  was  based  on  the  spin  study  data  (Ref.  4).  The  second  waB  sup¬ 
plied  'by  NASA  Langley  and  was  based  on  several  NASA  tunnel  investigations 
(e.g.,  Refs.  5  and  6).  The  Ref.  3  tabulation  was  a  function  of  a,  p,  and 
Batab  for  "the  razees 

0  <  a  <  90  deg  j  0  <  p  <  4o  deg  )  -21  <  &Bta,b  <  0  de® 

The  NASA  tabulation  was  a  function  of  a  and  p  for  the  range: 

-10  <  a  <  120  deg  j  -40  <  P  <  +40  deg 

Both  sets  are  limited  to  the  very  low  speed  regime  (M  =  0.2).  As  a  result 
of  comparison,  several  erroneous  data  points  were  discovered  in  the  Ref.  3 
data  and  were  reported  to  the  AFFDL. 

The  third  set  of  data  (Ref.  7)  was  obtained  from  the  NASA  Langley  Differ¬ 
ential  Maneuvering  Simulation  (DMS)  investigation  of  the  F-4j  and  slatted  F-4E 
aircraft.  Several  key  damping  derivatives  were  updated  from  Ref.  7  based  upon 
NASA  correlation  between  several  sets  of  tunnel  data  and  evaluation  of  aircraft 
response  obtained  in  the  simulation.  These  data  also  are  in  look-up  table  for¬ 
mat  with  coefficients  as  functions  of  a,  M,  and  H  over  the  range; 

-5  <  a  <  30  deg  }  .2  <  M  <  2.4  }  15,000  <  H  <  45,000  ft 

The  sideslip  coefficients  ere  valid  over  the  range  p  =  ±4o  deg.  Thrust  and 
drag  effects  are  also  modeled  in  detail. 
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As  originally  received,  the  data  did  not  coalesce  at  0.2  M,  The  data 
of  Ref.  4  were  therefore  given  the  most  weight  at  this  low  Mach  and  were 
smoothed  into  the  DMS  data  by  about  0.4  M.  The  resulting  data  provide  an 
excellent  model  of  the  aircraft  over  the  range: 

-4  <  a  <  30  deg  j  -10  <  p  <  10  deg  j  0.2  <  M  <  1 .0  ;  0  <  H  <  25,000  ft 

Besides  extension  of  coefficients  for  high  angles  of  attack  and  sideslip, 
several  changes  were  made  to  the  Ref.  l6  data  to  simplify  use  of  the  data. 
These  included  normalizing  the  data  for  trim  Mach  Number  (determined  from 
previous  digital  simulations  using  the  Ref.  l6  data  set),  using  a  single  alti¬ 
tude  (l 5,000  ft),  and  defining  values  for  the  second  set  of  the  parameters 
Cnp,  C^,  and  ACm. 

J.  AERODYNAMIC  COEFFICIENT  SATA 

The  aerodynamic  coefficients  are  tabulated  in  Table  9  on  pages  25  to  27 
and  plotted  in  Figs,  5  to  28.  All  coefficients  are  functions  either  of  a  or 
of  a  and  p.  The  aerodynamic  data  file  in  Table  9  includes  a  one-line  iden¬ 
tifier  for  each  coefficient,  plus  one  or  two  lines  describing  the  independent 
variable(s),  The  formats  of  these  identifying  lines  are  given  In  Table  8, 

All  coefficients  are  presented,  five  to  a  line,  in  5-degree  angle-of- attack 
increments . 

For  those  coefficients  which  are  functions  of  0  <  a  <  110°,  the  values 
given  correspond  to  a  as  follows: 

(0) . (20) 


(100)  (105)  (110) 

For  Oia<  45°  the  values  are  as  follows: 


(0) 

(25) 


(20) 

(45) 


TABLE  8. 

DATA  FILE  IDENTIFICATION  LINES 


Coefficient  identi¬ 
fier: 

VARIABLE 

NAME 

VARIABLE 

UNITS 

NUMBER  OF 
INDEPENDENT 
VARIABLES 

DATA  POINTS 
PER  LINE 

NUMBER 

OF  LINES 

TOTAL  NO- 
OF  POINTS 

(1  line) 

(Afl) 

(Afi) 

(12) 

(12) 

(12) 

(13) 

Independent  vari¬ 
able  identifier) a) i 

INDEPENDENT 

VARIABLE 

NAME 

INDEPENDENT 

VARIABLE 

UNITS 

MINIMUM 

VALUE 

INCREMENT 

MAXIMUM 

VALUE 

NO-  OF 
VALUES 

(1-2  linei) 

(AS) 

(AS) 

(B14-6) 

(E14- 6) 

(Ell*. 6) 

(12) 

The  coefficients  DCMl  and  DCM2  are  functions  of  both  a  and  p.  In 
this  case,  a  is  incremented  first,  then  p: 

(0,0) . (20,0) 

(25,0) . (1)5,0) 

(0,5) . (20,5) 


(25,30) . (a,  p) 

As  an  example  in  interpreting  the  data  file,  the  file  for  Cyp  Is 
given  below: 


CYB 

PER  DEG 

1  5  5  23 

ALPHA 

DEGREES 

0. ,5. ,110. ,23 

-.011 

-.011 

-.011  -.011 

-.011 

-.011 

-.011 

-.0092  -.0079 

-.0065 

-.0050 

-.0036 

-.0023  -.0008 

.0006 

.0020 

.0030 

.0030  .0030 

.0030 

.0030 

.0030 

.0030 

The  first  line  identifies  the  coefficient)  gives  the  units;  indicates 
that  it  is  a  function  of  one  variable,  and  that  it  is  listed  five  points 
per  line  in  five  lines,  with  a  total  of  ?3  data  points  given-  The  second 
line  identifies  a  as  the  independent  variable;  gives  its  units;  gives  the 
minimum  value  as  0  deg,  the  Act  increment  as  5  deg,  and  the  maximum  value 
as  110  deg,  for  a  total  of  25  separate  Act  points- 
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TABLE  9-  AERODYNAMIC  COEFFICIENTS 


S.O.  1081  EXTENDED  AO  A  HAT  A  »  CRDi  S  CRH2  CHANGED  11/20/78 

OLKAS  1  3  3  23 

ALPHA  DEGREES  0 .  » 3 .  »  1.10 .  ,23 

.1220  >4146  ,4834  .9053  .9396 

1.010  1.067  1.049  1.024  .949 

.92?  .872  .798  .477  .344 

.418  .239  .137  .023  -.114 

-.237  -.402  -.347 


CL8TAB 

PER  DEG 

1  5  5  23 

ALPHA 

DEGREES 

0.. 5. >110. *23 

.00499 

.00494 

.00474  . 00,637 

,00371 

.00313 

. 00303 

.00406  ,00443 

.00438 

.00414 

.00373 

■00371  . 00348 

.00324 

.00300 

.00300 

.00300  .00300 

,00300 

.00300 

.00300 

.00300 

CDBAS 

1  3  3  23 

ALPHA 

DEGREES 

0.  *5. *110. .23 

.033  .044  .130  .247  .382 
.310  .453  . /B4  .911  1,511 
1.123  1.245  1.332  1.403  1,451 

1.493  1.523  1.542  1.370  1,585  - 

1.400  1.583  1.347 

CY1)  PER  DEG  1  3  3  23 

ALPHA  DEGREES  0. *3. » 110. r23 

-.oil  -.on  -.on  -,<m  -.on 

-.011  -.011  -.0092  -.0079  -.0045 

-.0050  -.0034  -.0023  -.0009  .0004 

.0020  ,0030  .0030  .0030  .0030 

.0030  .0030  ,0030 

CYDR  PER  DEG  1  3  2  10 

ALPHAS  DEGREES  0. ,5**45. .10 

.00131  ,00142  .00134  ,00130  .00124 

,00098  .00073  .00048  .00024  0.0 

CRB1  PER  DEG  1  5  5  23 

ALPHA  DEGREES  0. ,5. >110. *23 

-.00132  -.00173  -.00234  -.00219  -.001.18 

-.00004  ,00023  -.00010  -.00108  -.0010 

-.00222  -.00241  -.00235  -.00248  -.00279 

-.00291  -.00300  -.00300  -.00300  -.00300 

-.00300  -.00300  -.00300 

CRB 2  PER  DEG  1  5  8  23 

ALPHA  DEGREES  0. *5 . r 110 . r 23 

-.00132  -.00173  -.00234  -.00219  -.00219 

-.00219  -.00219  -.00183  -.00152  -.00190 

-.00222  -.00241  -.00233  -.00248  -.0027V 

-.00291  -.00300  -.00300  -.00300  -.00300 

-.00300  -.00300  -.00300 


CRF'l 

PER  RAD 

4  3  2  10 

ALPHAS 

DEGREES 

0. *3. *43. *10 

-.30  4 

-.300  -. 

255  -.102  -, 

i  100 

-.230 

-.304  -. 

250  -.164  -, 

,100 

CRP2 

PER  RAD 

1  5  2  10 

ALPHAS 

DEGREES 

0, *5. *45. *10 

-.304 

-.300  -. 

260  -.240  - 

,  240 

-.240 

-.240  -. 

230  -.144  -, 

,  109 

CRR 

PER  RAD 

1  5  2  10 

TABLE  9 


(CONTINUED) 


ALPHA 2  DEGREES  0.*5.r4S.* 10 
.O'), 4  ,074  .108  .170  .318 
.2)37  .301  .200  ,00,4  .050 


CKO  A 

PER  DEG  1 

S  2  10 

ALPHAS 

DEGREES  0. 

*3. *45. *10 

.000448 

,000459 

.000431  ,000403 

.000379 

.000354 

,0003.1.2 

.0002:10  .000103 

.  000000 

CRDSP 

PER  DEG  1 

3  2  10 

ALPHAS  , 

DECREES  0. 

*3. *45, *10 

. 000 1 30 

.000140 

,  000093  .00004,4 

.000054 

.000040 

. 000030 

.000020  ,000010 

.000000 

CRtl  K 

PER  DEG  1 

5  2  10 

ALPHAS 

DEGREES  0. 

.•5,  *43.  r  10 

,000199 

.  0001413 

.000133  ,000114 

.000092 

.000049 

.000032 

.000033  ,000017 

.000000 

CM DAS  1  S  3  23 

ALPHA  DEGREES  0 .  *  3 .  *  1 .1 0 .  . 23 
.000  -.013  -.027  -.033  -.030 

-.0451  -.093  -.144  -.184  -.242 

-.2,43  -.21)3  -.310  -.381  -.430 

-.478  -  .301  -.332  -.371  -.414 

-.ASH  -.493  -.729 

DCM1  2  314  70 

ALPHAS  DEGREES  0 , » 5 .  *  4 S .*10 
BETA  DEGREES  0..3..3O.,  7 
.0000  .0000  ,0000  ,0000  .0000 

,0000  .0000  ,0000  .0000  .0000 


-.0000 

-,0080 

•■.0080 

- , 0138 

-.0099 

.0018 

,0005 

.0000  .0000  ,0000 

-.0170 

-.0170 

-.0170 

- . 0238 

- , 0506 

-.0204 

-,0160 

1000 

,  0000  . 

0000 

-.0312 

-.0312 

-.03151 

-.0480 

- , 0800 

•••.0443 

-.0420 

- .  0023 

.  0000 

,0000 

-.0433 

-.0435 

-.04158 

-.0705 

-.0800 

-.07,48 

- , 0755 

-.00(37 

.0000 

,0000 

-.0343 

- .  03  A  S 

-.03,45 

-.0080 

-.1000 

-.0933 

-.0944 

-.0106 

,0000 

.0000 

-.0474 

- .  0,474 

■ . 0674 

-.1071 

- . 1200 

-.1138 

-.1133 

-.0125 

.0000 

.0000 

DCM3  2  314  70 


ALPHAS 

DEGREES  0 , *  8 . 

*  45 , * 10 

BETA 

DEGREE 

8  0 .  ,  3 . 

*  30.  *7 

0,0  0, 

,  0  0.0 

0.0  0 

,0 

0.0  0 , 

,0  0.0 

0.0  0 

.0 

.  0080 

.  . 0080 

.0080 

.0138 

.0099 

-.0018 

■  -.0005 

.  0000 

.000 0 

.0000 

,0170 

.0:1,70 

.0170 

.  02  30 

,0306 

.0204 

.0160 

.0000 

,0000 

.0000 

.0312 

.0312 

.0312 

.0450 

,0580 

.04,43 

.04251 

.  0028 

,0000 

.0000 

.0435 

,0433 

» 0455 

,0705 

,0800 

.07,48 

.0753 

.  0087 

,0000 

.0000 

,  05,45 

.  0563 

,  0368 

.  0880 

,1000 

.  0933 

.0944 

.0106 

,0000 

.  0000 

,0,474 

.0674 

,0674 

.  107:1 

.1200 

.1138 

.1133 

.01213 

.0000 

.0000 

2.6 


TABLE  9 


(CONCLUDED) 


CMC! 

PER  RAD 

1  3 

2  10 

alpha:; 

DEGREES 

0 . ,  5 , 

,45.  ,10 

-3 , 10 

-3.41  -3 

,  78 

-3.92  -3 

.80 

-3.62 

-3 .  15  -3 

.13 

-3,115  -3 

.18 

CHAD 

PER  RAD 

1  B 

2  10 

ALPHAS 

DEGREES 

0  .  ,  5 . 

,45.  ,10 

-1,30 

-1.42  -1 

.  38 

-1.63  -1 

.59 

-1.52 

-1 . 32  -1 

.32 

-1.32  -1 

,32 

CMSTAB  PER  DEO  1  8  13  23 
ALPHA  DEGREES  0. *3 . , 110. ,23 
-.0:1010  -.01004  -.009131  -.00924. 

-.00743  - *00733  ".00677  -.00627 


-.00819  -.00 4 A A  -.00414  -.00361 

-.00256  -.00202  -.00130  -.00099 

0.00000  0.00000  0,00000 
GMDA  PER  DEG  1  5  2  10 
ALPHAS  DEGREES  0.. 3.  .45.  .10 


-.00103  -.00092  -.00079  -.00068 

- . 00066  -.00066  -.00066  -.00066 
CMDBP  PER  DEG  1  S  2  10 
ALPHAS  DEGREES  0,.S.# 48 .#10 


- .00828 
-.00572 
- . 00308 
- . 00032 


-.00066 

-.00066 


.000130  .000079  ,000033  .000007  .00000 
.000000  .000000  .000000  ,000000  .000000 
CM81  PER  DEG  1  8  2  10 
ALPHAS  DEGREES  0. #8, .48. #10 
. 0022  ,0022  .0021  ,0012  -.0006 

-.0033  -.0040  -.0038  -.0024  -.0020 

CNU 2  PER  DEG  1  S  2  10 

ALPHAS  DEGREES  0. . 3 . r  48 . . 10 
, 0022  .0022  .0021  .0012  .0000 

-.0003  -.000(3  -.0011  -.0014  -.0017 

LHP  PER  PAD  1  S  2  10 
ALPHAS  DEGREES  0. #3. .43. #10 
,003  -.014  ,007  .023  -.002 

. 088  .122  ,080  .041  .000 

CNR  PER  RAO  1  5  2  10 
ALPHAS  DEGREES  0,. 5. .43., 10 
-  .373  -.361  --.361  -.'370  •  ,302 

-.660  -.360  -.372  -.190  .000 

CHDA  PER  DEG  1  3  2  10 
ALPHAS  DEGREES  0 . #3 . » 48 . , 10 
-.000047  -.000040  000034  -.000067  -.000068 


-.000068  -.000069  -.000070  -.000070  -.000070 

CNDSP  PER  DEO  1  3  2  10 

ALPHAS  DEGREES  0 . , 5 . » 45 . » 10 

.0000332  .0000433  .0000280  ,0000132  .0000090 

,0000044  .0000009  .0000000  .0000000  .0000000 

CM DR  PER  DEG  1  3  2  10 

ALPHAS  DEGREES  0. #3. #43. *10 

-.00090  . .  -.00078  -.00073  -.00065 

-.00030  -.00032  -.00023  -.000:14  ,00000 
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Figure  5.  CLbasic  (a) 


Figure  6.  Cr  (a) 
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Figure  25.  Cn  (a) 


G.  VALIDATION 


No  F-4J  flight  test  traces  of  high  AQA  are  available  to  compare  against 
dynamics  from  the  foregoing  mathematical  model.  However,  the  model  was 
originally  developed  for  use  in  a  moving-base  piloted  simulation  .  •„  >'hich 
Navy  F-4J  pilots  were  given  air  combat  maneuvering  training.  As  a  part 
of  the  checkout  and  acceptance  tests,  the  simulation  was  flown  through  var¬ 
ious  offensive  and  defensive  combat  maneuvers,  stalls,  and  departures  by 
Navy  instructor  pilots  who  indicated  it  adequately  represented  the  F-4J 
handling  and  performance.  This  provided  the  first  gross  validation  of  the 
aero  model. 

The  only  high-AOA  flight  traces  available  are  from  the  F-4E  stall/  post- 
stall  flight  test  (Ref.  B).  McDonnell  Aircraft  Co,  indicates  that  all  hard- 
wing  models  (no  leading-edge  slats)  of  the  F-4  have  approximately  the  same 
stall/departure  characteristics.  The  Air  For^e  (Ref.  8)  noted  the  longer 
nose  F-J+E  to  have  somewhat  less  wing  rock  tendency  than  the  C  and  D  models. 
The  F-4j  is  more  similar  to  the  C  and  D  and  therefore  might  also  have  more 
pronounced  wing  rock.  Also,  the  F-4E  flight  test  vehicle  was  equipped  with 
a  spin  chute,  had  a  beefed-up  aft  fuselage  structure  to  handle  the  spin  chute 
loads,  and  had  offsetting  ballast  in  the  nose.  Thus  its  inertia  characteris¬ 
tics  are  quite  different  from  the  average  F-4J,  as  shown  in  the  following: 


F-4E  (Ref.  8) 

f-4j 

stores 

empty  tank 

empty  tank 

wing  stores 

pylons  1,2, 8, 9 

pylons  2,8 

W  (lb 

40,000 

37,000 

c.g.  (<f>  MAC) 

28.1 

29.3 

Ix  (slug- ft2) 

27,500 

23,850 

Iy  ( slug- ft2) 

157,000 

127,400 

Iz  (slug-ft2) 

180,600 

146,000 

Ixa  (slug- ft2) 

5,500 

2,  210 

This  difference  is  assumed  to  have  minor  influence  on  the  baGic  aerodynamic 
stall/separture  characteristics  as  a  function  of  AOA  but  a  significant 


42 


influence  on  dynamic  response  parameters  such  as  dutch  roll  frequency  and 
response  to  control  inputs.  With  this  in  mind,  we  will  proceed  with  a  com¬ 
parison  of  our  F-4j  model  against  the  F-4E  flight  test  results. 

Reference  8  indicates  that  the  clean  aircraft  exhibited  lateral/direc¬ 
tional  stability  breakdown  in  the  form  of  a  slightly  divergent  dutch  roll 
(wing  rock)  in  the  vicinity  of  19-22  deg  AOA.  As  AQA  was  further  increased 
the  motion  progressed  from  primarily  roll  to  yaw  (nose  slice)  in  the  region 
of  22-25.5  deg  AOA.  The  dutch  roll  mode  evaluated  over  the  same  AOA  region 
from  our  6  DOF  frozen  point  model  is  shown  in  Fig.  29*  This  also  indicates 
that  at  zero  sideslip  the  dutch  roll  slowly  becomes  divergent  at  about  19 
deg  AOA.  The  aero  data  plot,  Fig.  23,  shows  that  passes  through  zero  at 
about  20  deg  AOA,  while  at  this  same  point  (Fig.  10)  is  still  relatively 
large.  One  would  expect,  then,  that  the  dutch  roll  motion  would  be  primarily 
rolling  motion  (wing  rock) .  By  25  deg  AOA  C ^  is  very  small,  while  0 ^  is 
very  large  negative.  Thus,  one  would  expect  the  dutch  roll  mode  to  be  pri¬ 
marily  one  of  yaw  motion  (nose  slice). 

The  root  locus  of  Flg.'29  also  Indicates  that  the  dutch  roll  mode  is 
quite  sensitive  to  sideslip  in  this  same  AOA  range.  Therefore,  it  ip  neces¬ 
sary  also  to  look  at  time  trace  comparisons  between  flight  test  and  our  com¬ 
plete  6  DOF  nonlinear  aero  model.  Figure  30  is  an  F-4E  time  trace  from  Ref. 

8  at  21  deg  AOA  and  with  all  augmentation  off.  (it  also  has  tralllng-edge 
flaps  at  half  deflection  but  they  have  a  relatively  small  influence  on  lat¬ 
eral/directional  stability.)  The  traces  show  the  dutch  roll  to  be  a  constant- 
amplitude,  limit-cycle-like  oscillation.  A  similar  set  of  traces,  starting 
at  21  deg  AOA  from  our  F-4J  6  DOF  model  shows  an  almost  identical  oscillation 
(Fig.  31 ) .  A  comparison  of  specific  traces  shows  the  following  peak-to-peak 
excursions. 


Flight  Test 

Simulation 

a 

21  deg 

21  deg 

0 

1 1  deg 

6.3  deg 

<P 

30  deg 

21  deg 

P 

60  deg/sec 

24  deg/sec 

r 

2.5  deg/sec 

1 .4  deg/sec 

period 

3  sec 

6  sec 

maneuver 

wind-up- turn 

straight /level 
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Figure  29.  F-4J  6  DOF  Linearized.  Equation;  Lateral/ 
Longitudinal  Root  Migration  with  a  and  p 


Figure  30.  Controls  Fixed  Wing  Rock 
Time  History  (Ref,  8) 


The  absolute  magnitude  of  motions  depends  upon  the  initial  excitation 
and  is  not  particularly  close.  However,  the  ratios  of  the  motions  are  quite 
close  and  again  tend  to  validate  the  aero  model.  The  period  of  the  oscilla¬ 
tion  is  quite  different  because  of  the  higher  energy  (dynamic  pressure)  of 
the  flight  test  maneuver. 

The  nonlinear  response  traces  of  Fig.  3l  show  one  additional  factor. 

It  will  be  noted  that  the  longitudinal  traces  show  a  frequency  which  is  pre¬ 
cisely  double  the  dutch  roll  frequency,  while  Fig.  29  predicts  the  longitu¬ 
dinal  short  period  to  be  highly  damped  and  about  equal  in  frequency  to  the 
dutch  roll.  Other  F-4E  flight  test  time  traces  from  Ref.  8  (e.g.,  Fig.  32) 
reflect  the  same  "frequency  doubling"  noted  in  our  aero  model.  This  "fre¬ 
quency  doubling"  can  only  be  caused  by  the  pitching  moment  due  to  sideslip, 
Cmp,  which  is  quite  strong  and  negative  for  the  F-4  aircraft.  Note  in  Fig. 
32  that  each  peak  in  the  a  trace  coincides  with  a  zero  crossing  in  p.  In 
this  set  of  traces  the  rudder  activity  indicates  that  the  yaw  SAS  was  on, 

A  comparison  set  of  traces  from  the  F-4J  model  with  pitch  and  yaw  SAB  on, 
starting  at  a  trim  of  Oq  *  23  deg,  p0  «  o  is  shown  in  Fig.  33.  Again,  the 
motions  of  Figs.  32  and  33  are  remarkably  similar. 

A  final  comparison  between  flight  and  simulation  is  shown  in  Fig.  34. 
Both  aircraft  depart  from  a  wind-up  turn  to  the  left. 

On  the  basis  of  the  above  it  was  considered  that  the  aero  model  exhibits 
characteristics  adequately  representative  of  the  F-4  aircraft  for  use  in 
analysis  and  simulation  of  stall/departure. 
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Figure  34.  Departure  During  Left  Wind-Up  Turn 
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APPEND  DC  II 

1-1 4A  DATA  PACKAGE 


A.  OtTRODVOTION 

This  appendix  documents  the  mathematical  model  used  in  the  simulation 
of  the  P-14A  aircraft.  The  model  represents  the  "clean"  (no  slat  or  flap 
deflection)  aircraft  at  low  Mach  number.  All  aerodynamic  characteristics 
are  for  a  wing  sweep  of  22  deg  (Pig.  35)  which  'is  maintained  in  flight  by 
the  Mach-Sweep  Programmer  below  M  »  0.7.  Aerodynamic  coefficients  are 
derived  from  wind  tunnel  data,  with  modifications  based  on  parameter  iden¬ 
tification  using  flight  test  data. 

The  appendix  is  organized  as  follows:  the  equations  for  the  aerodyn¬ 
amic  forces  and  moments  are  presented  in  Section  B;  the  origins  of  the 
individual  aerodynamic  coefficients  is  discussed  in  Section  C;  the  coeffi¬ 
cient  "look-up"  tables  are  described  in  Section  D;  the  coefficient  data 
are  listed  in  Table  12  and  plotted  in  Figs.  %  through  71 . 

The  equations  of  motion  presented  in  Appendix  I,  Section  B,  are  applic¬ 
able  for  thiB  model, 

B.  AERODYNAMIC  FORCE  AND  MDMENT 

COEFFICIENT  EQUATIONS 

The  aerodynamic  model,  Table  10,  is  largely  based  on  an  F- 1 4  model  used 
on  the  NASA  Langley  Research  Center  Differential  Maneuvering  Simulator  (DMS), 
Ref.  9.  Modifications  were  made  to  incorporate  additional  wind  tunnel  data 
and  simplify  the  functional  dependence  on  sideslip  angle  where  appropriate. 
The  terms  ACl^so^  an^  were  added  to  allow  "tuning"  of  the  model 

using  flight-derived  trim  data.  Similar  berms  were  Included  in  the  lateral 
equations  to  allow  inclusion  of  nonzero  sidoforce,  yaw,  and  roil  moments  at 
p  a  0.  Such  moments,  arising  from  aircraft  asymmetries  and  asymmetric  vortex 
shedding,  have  been  suggested  as  contributors  to  departure  problems. 
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TABLE  10.  P-1 4  AERODYNAMIC  FORCE  AND  MOMENT  EQUATIONS 


UJ  A«rodyn»mlo  Force  Eguttlam 

01  “  ^BASIC*®'  P)  +  +  4Cl,gp(a)(  |Blp|/55)  +  AC^a) 


1  t&6§  <*>6« 

(8,  l 

>  -10°) 

1  -10CL8i/a^« 

+  (6,  ♦  10)C.  (a) 

(8,  < 

-10°) 

i(°)(  |p*p|/55>  + 

•i 

<6|  > 

-10°) 

“10cEfcg)(&)  + 

(6«  +  10)Cmt2(o) 

(6.  < 

-10°) 

‘V  ■  cyW8IC(<»»P)+  Cyp^j(a)  +  Cy8t(a)B»  +  Oy6r(o,  e)Br  +  cy6§p(“)6«p  +  (b/2Vo)[Cyr(o)r  +  Cyp(a)p] 
(b)  Aerodynwnio  Mdnwnt  Equation*  (Xog  In  %  g) 

0,11  ’  C“DA8IC(a'  p*  +  +  f<xc«-XMf)/l°o](Ci,  oat  a  +  Cj  iln  a)  +  AC^aK  |6ip|/35) 

*  +  (Sq/sVoJCnqfa), 

AC«nsg(a.)  «  (8,  >  -io°) 

AC^(a)  .  HOC^Ca)  +  (8,  ♦  10)0^,  (a)  (6,  <  -10°) 

°n  ■  C"BABIC(o'P)  +  +  ACn8|l(u)6«  *ln  (8,S  P)  +  t(Xeg  ~  W )/lOO](S/b)Oy  + 

+  Cn6r(o,p)6r  +  CnBip(o)6ip  +  <b/2V0)[c„r(a)r  ♦  Cnp(a)p] 

C 1  ‘  CiBA8IC^a'  P^  +  A%.o^  +  Cl6^a»  P)  »•  +  cl6j,(0»  P)  6r  ♦  Cl6ip(a^»P  +  (b/2V0)[C|r(o)r  ♦  CJp(o 

WE:  For  conciseness,  Bgtab  has  been  reduced  to  B 


b: 


l 
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la  the  DMS  model  the  static;  yawing  moment  coefficient  was  treated  as 
a  function  of  three  variables:  a,  p,  and  6a  (&s  =  &atab)*  The  unusual 
influence  of  stabilator  deflection  on  directional  stability  arises  from 
an  interaction  between  the  engine  inlet  shed  vortex  and  the  vertical  and 
horizontal  stabilizers,  Ref.  10.  Examination  of  the  data  indicated  that 
this  effect  could  be  modeled  as  a  yawing  moment  increment  which  was  a  func¬ 
tion  of  a,  Sa,  and  p  added  to  the  static  yawing  moment  Cn^siC .  Rather 
than  treat  this  increment  as  a  three  variable  look-up  table,  it  was  approx¬ 
imated  as 

ACnSa(a)5a  sin  (8.2p  ) 

where  A0n8  (a)  is  a  one  variable  look-up  table.  The  sin  (8.2p)  factor  in¬ 
sures  that  for  5a  <  0  (TEU)  the  effect  is  destabilizing  for  either  positive 
or  negative  sideslip.  The  longitudinal  stabilator  effectiveness  is  nonlinear 
in  Ba  in  the  DMS  model.  This  effect  was  retained  in  the  STI  model  but  the 
dependence  on  sideslip  was  eliminated. 

0.  DATA  SOURCES 

The  aerodynamic  data  package  employed  in  the  F-1^  simulator  was  put  to¬ 
gether  from  several  sources  and  then  modified  somewhat  to  obtain  an  accept¬ 
able  match  between  the  simulation  and  actual  flight  traces.  The  purpose  of 
this  section  is  to  document  the  originn  of  the  coefficients. 

Although  a  number  of  high  angle  of  attack  wind  tunnel  tests  have  been  con¬ 
ducted  for  the  F-1 4,  the  following  tests  were  the  primary  sources  for  the  STI 
model: 


15  March-l6  April  1971;  NASA  Ames  Research  Center, 

12  ft  pressure  tunnel 

August  1971;  NASA  Langley  Research  Center, 

30'  x  6o ’  (full  scale)  tunnel 

The  complete  test  reports  were  not  available  for  these  test  programs. 
Selected  ARC  data  were  evailable  in  Ref.  10  and  selected  LRC  data  in  Refs.  11, 
12,  and  13.  The  ARC  data  is  generally  considered  to  be  the  "best",  primarily 
because  the  Reynolds  number  is  higher;  however,  its  use  was  limited  in  that 
Ref.  10  contained  only  static  roll  and  yaw  coefficients  as  functions  of  a  and 
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p  and  normal  force  and  pitching  moments  as  a  function  of  a  only.  No  dynamic 
derivative  or  control  effectiveness  data  is  available  from  the  ARC  test. 

These  were  obtained  from  the  Ref.  9  DM3  simulation  which  was  based  on  the 
liRC  30'  x  6o'  data  with  some  modifications  based  on  simulation  pilots'  opin¬ 
ions.  The  DMS  data  appears  to  be  primarily  "raw"  LRC  wind  tunnel  data  as  it 
shows  characteristic  scatter,  and  coefficients  (such  as  Cn(a,p))  which  are 
not  symmetrical  in  p.  When  this  data  was  used,  coefficient  values  at  posi- 
tive  and  negative  p's  were  averaged,  where  appropriate,  to  produce  tables 
symmetrical  in  p.  In  some  cases  the  averaged  curves  were  further  smoothed 
to  reduce  questionable  scatter. 

The  adequacy  of  the  wind  tunnel  derived  aerodynamic  data  was  validated 
by  comparison  of  6  DOF  analytic  models  with  available  F-'l 4  flight  test  data. 

On  the  basis  of  these  comparisons  several  aerodynamic  coefficients  were  fur¬ 
ther  modified  as  explained  in  Section  E. 

Plots  of  the  coefficients  are  shown  in  Figs.  36  to  71  •  The  reference 

point  for  all  moment  coefficients  is  on  the  FRL  at  l6#  MAC.  The  spoiler 

produces  no  effect  above  a  =»  10  deg  and,  since  primary  interest  is  in  the 
high  a  regime,  all  spoiler  functions  are  zero. 

D.  AERODYNAMIC  COEFFICIENT  DATA 

For  use  by  the  simulation  program  all  aerodynamic  coefficients  are  repre¬ 
sented  as  piece-wi3e  linear  approximations  in  "look-up"  tables.  These  tables 
are  functions  either  of  a  only  or  of  a  and  p  (see  Table  11 ) .  In  either  case, 

the  a  breakpoints  are  every  5  deg  from  a  =  0  to  55  deg.  The  p  breakpoints 

are  every  5  deg  from  p  =  0  to  ±20  deg.  A  listing  of  the  computer  data  file 
containing  the  lookup  tables  is  shown  in  Table  12.  The  order  of  the  individ¬ 
ual  coefficient  tables  in  the  aerodynamic  data  file  is  as  shown  in  Table  11. 
For  functions  of  a  only,  a  is  incremented  by  column  and  then  row,  as  indicated 
below: 

0,  5  deg,  10  deg,  15  deg,  20  deg 
25  deg,  30  deg,  35  deg,  4o  deg,  45  deg 
50  deg,  55  deg,  0,  0,  0 


For  functions  of  a  and  0,  0  is  incremented  first,  then  a,  as  follows: 

Row 

1  (0,  -20),  (0,  -15),  (0,  -10),  (0,  -5),  (0,  0) 

2  (0  ,  5),  (0,  .10),  (0,  15),  (0,  20),  (5,  -20) 

22  (55,  10),  (55,  15),  (55,  20),  (o),  (o) 

For  convenience  in  creating  computer  mnemonics,  the  control  variable 
8  +  .  has  been  abbreviated  as  8a  in  the  data  file  and  plots. 


TABLE  1 1 .  ORGANIZATION  OF  AERODYNAMIC  DATA  FIDE 


PARAMETER 

function  or 

PARAMETER 

UNITS 

BTI  MNEMONIC 
AND  ARRAY 
dimension 

NUMBER  OF 
MTA  PTS 

CARD 

(USE) 

CODE 

CLBASIC 

*>  B 

— 

CIBAS(12,9) 

108 

BTI 100 

ACi^mq 

a 

- 

DCLBO( 12) 

12 

101 

°l6,. 

a 

ciDSi(ia) 

18 

10U 

°Kia 

a 

t/d*B 

CDD8S( ia) 

12 

105 

aci,»? 

a 

— 

DCLSP( 12) 

12 

106 

°HSASIC 

BE,  B 

— 

casAS(ia,9) 

108 

800 

°D»., 

a 

I/4** 

00081(12) 

IS 

SOU 

co»,a 

a 

1/lU 

cdc*2(i8) 

18 

805 

A0d» 

a 

— 

DCDSP(ia) 

18 

ao6 

Atbasic 

a,  B 

— 

niAstia,9) 

108 

1  >00 

ACrB-O 

a 

— 

uaao(ia) 

12 

301 

% 

a 

lA*d 

«K  ia) 

18 

503 

cVr 

a 

1/pa4 

ends) 

12 

304 

CP6a 

a 

1/d«« 

18 

303 

Cy»P 

a 

l/d«« 

01082(12) 

12 

306 

°»»r 

ACig.Q 


««(ia,9) 

ORBAS( 12,9) 
00880(18) 


OiP  a 

l/r«d 

CTO<18) 

12 

405 

Ojr  a 

l/wd 

onR(ia) 

12 

4o4 

0|6»  »*  9 

l/d«« 

aRnA(i8,9) 

108 

405 

°1B8P  * 

l/d«* 

CRDSP(IS) 

12 

406 

C«S-  °*  P 

l/d»« 

CRM(18,9) 

100 

4o? 

CMBAB(12,9) 

108 

300 

c“«ASIC  ’  p 

ACm.o  a 

— 

D0M80(  12) 

12 

501 

feq  » 

1/rtd 

0M9(18) 

12 

303 

l/d*« 

CMD81  ( IS) 

18 

304 

0*6, 8  a 

i/d  H 

CMDS2(12) 

12 

305 

MB8j>  a 

DCMaP(IB) 

12 

>06 

0nBA8IC  *•  9 

— 

ums(i8,9) 

108 

600 

« 

— 

DOHBO(ia) 

12 

601 

^HB,  ° 

l/d«* 

DCHD8(  12) 

12 

602 

Onp  a 

l/r»d 

CTfR(  12) 

12 

6o> 

Onr  a 

l/r*d 

CNH(  18) 

12 

604 

°"6» 

1/di* 

CWDA(  12) 

12 

603 

C>1BSP  “ 

l/d*i 

CNDSP(l2) 

12 

6o6 

On»r  1  a*  9 

l/d»* 

CN»(I2,9) 

106 

607 

.  ,M,;i  viutt-i  v  -v 


TABLE  12.  AERODYNAMIC  DATA  FILE 


O.4145O0OE-C1 
0. 14SS500E4-C0 
0.5744194E4-00 
0 . 4024052E4-00 
0. 10017S4E4-01 
0  *  980Q749E4-00 
0. 1274219E4-01 
0 i 121 3586E4-01 
0  *  1403140E4-01 
0 » 14404S0E4-01 
0 • 17442B9E4-01 
0 . 1790444E4-01 
0 . 1832S92E4-01 
0. 1840752E401 
0. 1770324E4-01 
0. 1748122E4-01 
0. 1464409E4-Q1 
0.  1448984E4-01 
0.1474011E4-‘01 
0. 1546793E+01 
0.13Q5247E4-01 
0.1394S44E4-01 
0.,0.»0.,0.»0., 
0.,0.,Qt>0.,0,, 


O.147OSO0E4-OO 
0.109B300E4-00 
0.4024032E4-00 
0.3744194E4-00 
0.1O19434E4-O1 
0.9215549E4-00 
0.1278946E4-Q1 
0.1439421E+01 
0. 1549673E4-01 
0  * 1498947E4-01 
0.1740003E+01 
0.1052892E4-01 
0.1790444E4-01 
0  > 1870520E+01 
0  •  1733031E4-01 
0 . 1771S12E4-01 
0.1347197E4-01 
0,  US2013E+01 
0 « 1499900E4-O1 
0.  1S21319E4-01 
0.1394564E+01 
0.138324784-01 


Q. ,0. ,0. >0. >0. r _ 

.0145, .0123, .0135, ,0150, .0163, 
,01.45,  .0130,  .0113,  ,0043,  .0040 » 

.0070, ,0043,0. ,0» ,0.  _  _ 

,  00 A3,  ,  0095 , .  0 i  It)  77oF2S7»  01 20 , 
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E.  VALIDATION 


This  section  documents  the  validation  of  the  F-l4  high  angle-of- attack 
aero  model  against  flight  test  data.  Basically,  it  was  necessary  to  deter¬ 
mine  whether  the  simulation  model  would  produce  the  high  a  dynamic  phenomena 
of  the  unslatted  F-14A,  specifically: 

•  Mildly  divergent  dutch  roll  in  the  AOA  range 
12  to  22  deg 

•  Wing  rock  when  AOA  is  maintained  in  the  12  to 
2 5  deg  range  at  low  Mach  Number 

•  Roll  reversal  starting  at  a  =  l8  deg  when  dif¬ 
ferential  tail  is  used  to  roll  the  aircraft 

1 .  Validation  Methods 

Two  basic  methods  were  available  for  comparing  the  simulation  model  with 
the  flight  test  data.  First,  simulation  time  response  traces  were  compared 
directly  to  the  flight  test  traces.  This  approach  had  limited  success,  since 
the  STI  unpiloted  time-domain  simulation  program  does  not  have  provision  for 
non-zero  acceleration  initial  conditions  (i.e.,  the  aircraft  is  assumed  to 
start  from  a  steady  trim  condition)  and  does  not  have  the  capability  to  re¬ 
produce  the  complex  control  inputs  of  the  flight  traces. 

The  second  method  involved  extracting  values  of  damping  ratio,  £,  and  fre¬ 
quency,  (u,  by  locally  approximating  oscillations  in  the  flight  traces  with  a 
second-order  linear-system  response.  In  regions  where  the  control  inputs  were 
fixed  or  could  be  related  to  identifiable  feedback  loops,  flight- derived  £ 
and  co  values  could  be  compared  to  open-or  closed-loop  roots  obtained  from  small- 
perturbation  linearization  of  the  simulation  model.  This  second  method  was 
particularly  useful,  since  the  lateral  oscillatory  response  was  of  primary 
interest  and  the  linear  response  parameters  £  and  co  can  be  directly  related 
to  the  aircraft's  stability  and  control  derivatives  through  literal  approxi¬ 
mate  factors.  Ref.  1 .  The  approximate  factors  were  rewritten  for  body- center¬ 
line-axis  stability  derivatives,  Table  13,  rather  than  the  usual  stability- axis 
factors,  so  that  the  response  parameters  could  be  related  directly  to  the  aero- 


dynamic  model  body  axis  moment  coefficients.  Thus  the  significant  stability 

and  control  derivatives  for  a  given  response  parameter  could  be  determined 
and  modified  to  match  the  flight  test  result.  Such  parameter  modifications 
made  by  linear  analysis  were  then  checked  by  extracting  £  and  o>  values  from 
non-linear  simulation  time  responses. 


TABLE  13.  APPROXIMATE  FACTORS  FOR  KEY  LATERAL - 
DIRECTIONAL  PARAMETERS  AT  HIGH  AOA  USING  BODY 
CENTERLINE  AXIS  DERIVATIVES  (y0  =  O,po  "  0) 


2  1  1 

=  Np  -  Lp  tan  8^ 


'  /  '  g  \  1 

2U(41  =  “Yv  “Lp  +  INp  “  tan~0o 


Np  -  Lp  tan  G0 
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i  L  -  -a 

%  \  VT, 


+  Nr  tan  8C 
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%  4  J|\|;{w-Van  e°>  -Np(4-Li  ta" 


2.  Flight  Test  Conditions  and  Configurations 

F- 1 L-  flight  test  data  covering  the  Wgh-a  region  was  obtained  from  the 
manufacturer,  Ref.  15.  From  thiE  data  six  flight  test  runs  (Table  1*0  were 
found  which  were  suitable  for  comparison  with  the  STI  simulation.  Three  dif¬ 
ferent  aircraft  were  involved  in  these  tests.  One  flight,  No.  199,  was  only 
used  for  comparison  with  the  simulation  trim  routine. 

Reproductions  of  the  flight  test  traces  are  shown  in  Figs,  7^,  75,  78, 

80,  82,  and  8^.  Aircraft  2  and  3  had  flight-test  nose  booms  with  angle-of- 
attack  and  sideslip-angle  sensors.  The  a  traces  for  Flights  230-5  and  230-6 
(Aircraft  3)  are  "raw"  noseboom  a  (denoted  ajjgD  on  the  traces).  Noseboom  a 
is  related  to  true  a  by 

^irue  =  *8561  Noseboom  +  •1'ir0  des 


82 
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Aircraft  IX  had  no  nose  boom  and  hence  no  sideslip  sensor.  Angle  of 
attack  was  measured  by  the  ARI  nose  probe  and  denoted  on  the  traces. 

The  nose-probe  a  is  related  to  true  a  by 

aTRUE  "  *8 122  aNOSEPRQBE  +  *797  deg 

The  simulation  model  was  designed  to  represent  the  "clean"  (flaps 
and  slats  retracted)  P~l4  at  low  Mach  number  and  forward  wing  sweep  (22 
deg),  Pig.  35.  The  flight  test  runs  used  for  validation  match  these  con¬ 
ditions  as  closely  as  possible,  specifically: 

•  Flaps  and  slats  are  retracted  except  for 
Bflap  =  —2  deg  on  Plight  No.  236-1 

•  Mach  number  below  M  =  .35 

•  Wing  sweep,  A  =  1(>“20  deg 

•  Glove  vane  and  speed  brakes  retracted 

•  No  external  stores 

The  maneuvers  were  performed  at  nz  0  1  g,  either  by  maintaining  speed  and 
AOA  constant  and  exciting  lateral  motions  with  a  doublet,  or  in  a  "1  g  stall" 
In  which  angle  of  attack  is  steadily  increased  as  speed  is  decreased.  Longi¬ 
tudinal  acceleration,  Vj,  varied  from  0  to  -0.4  g,  and  large  negative  flight 
path  angles  were  reached  at  higher  angles  of  attack.  Fairly  large  excursions 
in  mean  bank  angle  occurred  in  some  runs,  but  lateral  accelerations  were  gen¬ 
erally  small.  True  AOA  ranged  from  13  deg  to  32  deg,  which  covers  the  range 
of  interest  for  high  AOA  phenomena  on  this  aircraft. 

At  low  speeds  the  lateral  SAS  gains  are  scheduled  with  AOA  (Figs.  72  and 
73)  such  that: 

•  yaw  SAS  off  a  >  10° 

•  Blatst,ick  Br  Phasod  in  10°  <  a  <  20° 

•  Roll  SAS  and  command  20°  <  a  <  3^0 

augmentation  phased  out 


04 


N 


Stabler 

Actuator 


Rudder 

Actuator 


Where  a  as  used  in  Pigs.  72  and  73  is  the  angle  of  attack  as  measured  by  the 
ARI  nose  probe.  According  to  this  schedule  the  yaw  SAS  should  be  off  for  all 
cases.  However,  Flight  No.  243  appears  to  have  the  yaw  SAS  on  as  will  be  dis¬ 
cussed  later.  Plights  with  roll  SAS  and  ARI  on  were  included. 

3.  Flight  Teat  Comparisons  and  Derivative 
Modifications  to  Match  "Wing  Rook" 

Flight  199  (Aircraft  No.  2) 

Flight  199 >  Pig.  74,  is  a  1-g  stall  maneuver  in  which  angle  of  attack  is 
increased  steadily  from  10  deg  to  almost  40  deg  with  a  ramp  (TEU)  horizontal 
stabilizer  input.  This  maneuver  was  used  to  check  the  simulation  trim  values. 

At  a  given  speed,  1~g  trim  values  of  a  and  8atab  from  the  simulation  were  com¬ 
pared  to  the'  flight  a  and  8stab  traces.  Initial  comparison  using  the  original 
wind  tunnel  derived  aerodynamic  model  gave  good  comparison  with  flight  a  but 
the  flight  bBtab  traces  indicated  a  more  trailing-edge-up  deflection.  The 
original  simulation  Cra  (a,p)  was  based  on  the  LRC  wind  tunnel  data  (which  dif¬ 
fered  from  the  ARC  data  by  an  essentially  constant  nose  up  increment) .  A 
ACm^Q^)  increment  of  -.055  was  added  to  the  simulation  to  bring  the  pitching 
moment  into  line  with  the  ARC  data.  This  change  resulted  in  a  better  match  of 
the  bstab  trace,  as  shown  by  the  circled  points  in  Fig.  74. 

Flight  243  (Aircraft  No.  3) 

In  the  time  "slice"  of  interest  (t  «  33-50  sec)  AOA  is  maintained  at  l8  ± 

2  deg,  Fig.  75.  In  this  AOA  region  the  gain  schedule  implies  that  the  yaw  SAB 
and  ARI  are  off  and  that  the  roll  SAS  is  on.  However,  comparison  of  the  yaw 
rate  and  rudder  position  traces  indicates  that  the  rudder  is  correlated  with 
yaw  rate  at  approximately  the  magnitude  anb  phase  angle  that  would  result  from 
the  yaw  damper.  The  rudder  pedal  and  lateral  stick  are  both  effectively  zeroed, 
and  thus  the  aircraft  response  can  be  considered  the  free  response  of  the  aug¬ 
mented  airframe . 
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Figure  74.  Flight  199  Traces 
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The  aircraft  exhibits  a  divergent  "wing  rock"  oscillation.  The  fact 
that  the  envelope  appears  to  grow  linearly  with  time  rather  than  exponentially 
indicates  a  somewhat  nonlinear  response;  however,  approximate  values  of  damp¬ 
ing  ratio  and  frequency  can  be  extracted  as: 

t  =  -.05 

a  =  1 .1+  rad/sec 

The  dutch  roll  root  obtained  from  linear  analysis  using  the  original 
(wind  tunnel  derived)  aerodynamic  data  is  stable  and  about  15  percent  lower 
in  frequency  than  the  flight  value.  This  is  shown  (both  SAS  off  and  SAS  on) 
in  Pig.  76*  The  body-axis  approximate  factors,  Table  1 3 ,  evaluated  with  the 
original  aerodynamic  data  yield 


(-  tan  e0  Lp  j 1  '/£' 

• 

C 

(-  ‘“SoTr 

■ 

1.21  rad/sec 

^d 

t[”yv  “  lp  +  twre0-  (nj.  -  vf^)j 

ss 

4  (.0^75  +  .5758  -  .097O  -  ,5510) 

* 

.0863 

These  approximations 

are 

within  b  percent  of  the  values  obtained  by  numer. 

ical  factorization  of  the  characteristic  equation.  It  can  be  seen  that 
either  an  increase  in  dihedral  effect,  or  a  decrease  in  roll  moment 

of  inertia,  Ix,  would  increase  the  dutch  roll  frequency.  A  reduction  in 
the  magnitude  of  the  roll  damping,  C£p,  or  a  more  positive  value  of  the 
cross  coupling  parameter,  Cnp,  would  reduce  the  dutch  roll  damping.  Such 
changes  were  made  iteratively  using  the  approximate  factors  for  guidance 
until  the  dutch  roll  reasonably  approximated  the  observed  frequency  and 
divergence  rate. 
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Figure  T^,  Comparison  of  Oscillation  Root  from 
Simulation  v/lth  Flight  Tost  for  Flight  No.  245 
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These  changes  produced  a  small-perturbation  (linear',  dutch  roll  root  at 
Kd>  =  ["  *071;  1.4l]  for  the  bare  airframe,  and  [  £a,  tu^]  =  [— .04l, 

1.43],  BAS  on,  which  can  bn  compared  with  the  flight-derived  root  in  Fig.  76. 

To  check  possible  nonlinear  effects  the  airframe  alone  "Igital  simulation 
was  run,  Fig.  77,  with  the  final  modified  parameters.  The  pproximate  £  and 
in  in  response  to  s  pulse  rudder  were  extracted  as  for  the  flight  test  data. 
This  resulted  in  [£,u>]  =  [-.045,  1  *4]  for  the  bare  airframe  Actually  the 

frequency  started  at  about  1.45  rad/sec  and  decreased  to  about  1.35  rad/sec 
as  the  amplitude  of  oscillation  increased.  With  the  SAS  on  the  root  moves  to 
[£,w]  =  [-.024,  1  .40]. 


The  Flight  234-1  maneuver  was  performed  at  an  approximately  constant  AOA, 
Oq  =  15  deg,  Fig.  78*  In  the  time  slice  of  interest,  t  =  15  to  28  seconds, 
the  yaw  SAS  is  off  per  the  a  schedule,  as  can  be  confirmed  from  the  rudder 
position  trace.  The  ARI  is  scheduled  at  60  percent  but  has  no  effect  since 
the  lateral  stick  is  fixed.  The  roll  SAS  which  would  normally  be  on  at  this 
AoA  has  been  switched  off,  as  can  be  seen  from  the  roll  series  actuator  trace. 
Thus  from  t  =  1 9  to  L3  seconds  the  response  is  that  of  the  bare  airframe. 

This  shows  an  oscillating  divergence  (dutch  roll)  with  [£,^]  =  [-.07,  1.4  rad/ 
sec] . 

Calculation  of  the  dutch  roll  root  using  the  original  (wind  tunnel)  aero¬ 
dynamic  data  again  produced  a  stable  mode.  Estimates  of  parameter  changes 
required  to  match  the  flight  test  divergence  indicated  that  a  level  cf  roll 
damping  comparable  to  the  final  Flight  243  match,  C»  =  -.08  1 /rad,  was  re- 
quired.  When  this  C_gp  value  was  used  in  the  nonlinear  simulation  of  the  bare 
airframe.  Fig.  79,  the  approximate  damping  ratio  and  frequency  were  found  to 
be  [£,o>]  =  [-.058,  1.43]. 


The  Flight  236-1  maneuver  was  a  1  g  stall  in  which  angle  of  attack  was 
steadily  increased  from  a  =  13  to  28  deg  in  the  time  slice  examined.  Fig.  00. 
The  roll  SAS  has  been  switched  off,  as  seen  from  the  roll  series  actuator 
trace  and  the  lateral  stick  is  fixed,  thus  the  differential  tail  is  zeroed. 

The  yaw  SAS  is  also  off  per  the  a  schedule.  Between  t  =  35  and  40  sec  the 
lateral  stick  and  rudder  are  used  to  excite  a  rolling  oscillation.  The 
lateral  stick  is  then  centered  and  held  fixed.  A  rudder  pedal  input  is  made 
at  t  =  43  sec.  Also  note  that  the  rudder  apparently  was  used  to  counter  some 
other  disturbance  since  the  aircraft  is  rolling  to  the  right  and  the  rudder 
input  is  trailing-edge  left.  Thus  the  rudder  input  and  other  disturbance  may 
contaminate  the  response  sufficiently  to  alter  the  apparent  £  and  o>  values. 

As  in  the  Flight  243  trace,  the  envelope  of  the  bank-angle  oscillation  appears 
to  grow  linearly  with  time  rather  than  logarithmically,  but  an  estimate  of  the 
response  parameters  gives  [£,<n]  =  [-.085,  1  .30  rad/sec] . 

Matching  of  this  flight  condition  was  done  directly  from  nonlinear  simu¬ 
lation  time  responses,  Fig.  81 ,  Matching  frequency  required  a  30  percent  in¬ 
crease  in  lx  (at  this  heavy  weight  the  aircraft  has  approximately  4,000  lb  of 
fuel  in  the  wings)  and  matching  the  divergence  rate  required  reducing  the  roll 
damping  below  the  level  of  the  wind-tunnel  data  in  the  20  <  a  30  deg  region 
(again  similar  to  the  0 ^  reduction  required  at  15  <  a  <  20  deg).  The  final 
simulation  time  response  shows  a  nonlogarithmic  q>  envelope,  similar  to  the 
flight  traces,  which  is  possibly  a  nonlinear  effect  of  the  time-varying  a. 
Oscillation  frequency  derived  from  the  simulation  cp  trace  varies  from  1 .40  rad/ 
sec  at  low  amplitude  to  1.1  rad/sec  at  large  amplitude.  Damping  ratio  varies 
from  -.1 1  to  -.085. 


Flight  230-* 


In  the  time  slice  examined,  the  angle-of-attack  was  roughly  constant  at 
a  =  29  deg.  Fig.  82.  Although  there  is  considerable  stick  activity,  the  roll 
SAB  and  command  augmentation  is  almost  phased  out  at  this  AOA  and  thus  the 
differential  stabilizer  is  essentially  zeroed.  The  large- amplitude  rudder 
oscillation  Is  due  primarily  to  the  lateral  stick  deflection  acting  through 
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Figure  80.  Flight  256-1  U 
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Figure  8l .  Bare  Airframe  Simulation  Time 
Responses  for  Flight  236-1 


the  atick-to-rud&er  interconnect  (SRI).  From  an  examination  of  the  rudder 
and  roll  traces  it  cannot  be  determined  immediately  if  the  pilot  was  attempt¬ 
ing  to  excite  the  roll  with  the  rudder  or  was  feeding  p,  cp,  or  both  to  rudder 
to  damp  the  oscillation.  The  latter  possibility  can  be  checked  by  consider¬ 
ing  a  system  survey  (root  locus  and  Bode-siggiej  see  Part  I,  Section  II)  for 
<p  6r  loop  closure  as  shown  in  Fig.  83.  It  can  be  seen  that  at  gains  on 
the  order  of  those  observed  in  the  flight  traces,  =  1  deg/deg,  the  dutch 
roll  is  stable  but  a  divergent  first-order  root  appears.  Since  no  aperiodic 
divergence  is  apparent  in  the  flight  traces,  it  would  seem  likely  that  the 
rudder  input  is  an  open-loop  excitation.  From  the  flight  traces,  the  fre¬ 
quency  is  to  =  1.4-1 .6  rad/sec  and  |<p/8r  (ju>)|  =  .6  to  1.0  (-4.4  dB  to  0  dB). 
This  is  in  reasonable  agreement  with  the  open-loop  Bode  plot  of  Fig.  83  indi¬ 
cating  that  the  rudder  effectiveness  is  reasonable  at  this  angle  of  attack, 
a  =  29  deg. 

Flight  230-6 

From  t  =  20  to  4o  seconds,  AOA  is  increased  rapidly  with  a  ramp  stabila- 
tor  input,  Fig.  84,  The  rudder  and  lateral  stick  are  essentially  zero  during 
most  of  this  time.  There  is  no  evidence  of  a  divergent  lateral  oscillation, 
which  is  in  general  agreement  with  the  simulation  which  indicates  a  stable 
airframe  in  this  angle-of-attack  region  (a  >  25  deg).  The  simulation  yields 


a 

25 

+  .097 

1 .22  rad/sec 

30 

+  .244 

.93  rad/sec 

4.  Roll  Reveraal  Validation 

Figure  85  shows  the  migration  of  the  u*p  root  with  AOA.  It  can  be  seen 
that  between  a  =  15  and  20  deg  sty  becomes  real  with  one  zero  in  the  right 
half-plane.  This  produces  an  initial  open-loop  roll  opposite  to  that  com¬ 
manded  with  lateral  stick.  It  also  results  in  a  first-order  lateral-direc¬ 
tional  divergence  if  the  cp  —*■  loop  is  closed,  Fig.  86.  This  behavior  can 
be  traced  to  the  change  from  proverse  to  adverse  yaw  at  a  =  18  deg  shown  in 
the  Cng  vs.  a  plot  of  Fig.  87. 


(a),  l/deg 
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